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ABSTRACT 

We report on time-dependent axisymmetric simulations of an X-ray excited 
flow from a parsec-scale, rotating, cold torus around an active galactic nucleus. 
Our simulations account for radiative heating and cooling and radiation pressure 
force. The simulations follow the development of a broad bi-conical outflow 
induced mainly by X-ray heating. We compute synthetic spectra predicted by 
our simulations. The wind characteristics and the spectra support the hypothesis 
that a rotationally supported torus can serve as the source of a wind which is 
responsible for the warm absorber gas observed in the X-ray spectra of many 
Seyfert galaxies. 
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1. Introduction 

Approximately 50% of low-red-shift AGNs exhibit spectra rich in atomic absorption 
features in the X-ray band. These spectra, referred to as warm absorbers, contain numerous 
lines and bound-free edges from ions of intermediate- Z elements, and indicate blue-shifts 
of the absorbing material of ~ 10^ — 10^ km s~^. The term refers to the fact that the 
X-ray absorbing gas has an electron temperature and ionization degree which are interme- 
diate between that expected for neutral and fully ionized material. Warm absorber gas has 
been known for more than two decades, beginning with observations using the Einstein Ob- 



servatory satellite (Halpern 1984), although the limited spectral resolution of early X-ray 



experiments hampered a detailed study of its properties. 
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Recent X-ray telescopes, Chandra and XMM-Newton, provide unprecedented spectral 
resolution up to ~10 keV, and have revealed much of the complexity of warm absorber 
spectra. The ubiquity of X-ray warm absorbers and the implication of outflow in objects 
which are likely powered by accretion, combine to make warm absorbers of great interest for 
the understanding of active galactic nuclei. 

While observationally the reality of the X-ray absorbing gas is firmly established, the 
source of this gas remains uncertain. If the wind originates from an accretion disk 
(Konigl &: Kartje 1994 Murray et al. 1995), then warm absorber gas can occupy a vast spa- 



tial region from the inner accretion disk up to and beyond the broad line region. In addition, 
it is likely to have a speed which is comparable to the escape velocity in the region where the 
wind originates. If so, the line widths favor a wind originating far outside the inner accretion 
disk, at ~ 1 pc, for an AGN containing a 10^ black hole. 

The natural source of gas at this distance is the molecular torus which is responsible 
for obscuring the broad line region in Seyfert 2 galaxies, and which is thought to exist 
in some form in many low and intermediate luminosity AGN ( lAntonucci fc Miller 1986|) . 



The existence of an outflow from the torus has been suggested by (Krolik &: Begelman 1986 



Krolik fc Begelman 1988), and as the source of warm absorber flows by flKrolik fc Kriss 1995 



IKrolik fc Kriss 200ip . 



Direct interferometric observations of this dusty structure have recently become avail- 
able: mid-infrared high spatial resolution studies of the active nucleus of NGC 1068, revealed 
a dust torus - like structure that is 2.1 pc thick and 3.4 pc in diameter (IJaffe et al. 2004]) . 



Evidence of at least two temperature components has been established: the flrst one is ex- 
tended and warm (320 K) and conceals the second, which is compact and hot (> 800K). It 
is natural to assume that the latter component can be attributed to the inner part of the 
torus, heated by the radiation from the compact nucleus. Although such evidence is most 
solid for nearby active galaxies it is likely that the same obscuring torus paradigm may ap- 
ply in quasars whose central regions are heavily obscured by gas and dust (Type II quasars) 
(IZakamska et al. 2006p . 



The work done so far in modeling torus warm absorber flows employ either one-dimensional 
models e.g. (IChelouche fc Netzer 2005]) or two dimensional models which omit rotational 
forces and also employ an artificial shape for the torus, which is adopted as a boundary 
for the flow (IBalsara fc Krolik 1992p . These calculations support the idea that the outflow 
originates at r ~ 1 pc - the estimate supported both from virial arguments and also deduced 
from the warm absorber line widths. The content of the warm absorber gas is highly ionized 
cosmically abundant heavy elements. This gas can occupy a region from 0.01 pc to ~ 10 pc 
from the nucleus. The column densities are in the range 10^^ — lO^^cm"^. These estimates 
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are consistent with constraints on the position of the warm absorber gas, which are deduced 
from the absence of correlated variabihty of the hne equivalent widths to changes of the 
continuum (INetzter et al. 2003^ iBehar et al. 20031) . 



The goal of the current work is to show that a state-of-the-art hydrodynamic model 
of the evaporative torus wind possesses key characteristics of the the warm absorber. To 
approach this goal, in this paper we present results of models containing the following in- 
gredients: 1) Time-dependent axisymmentric numerical simulations that take into account 
rotational forces , X-ray heating and cooling and radiation pressure both in spectral lines 
and continuum; 2) Exploration of the resulting distributions of density, velocity etc.; and 
3) A calculation of the synthetic X-ray absorption spectra using the density and velocity 
distributions as input. Although calculating the spectra is crucial to check the relevance of 
the particular hydrodynamic solution to the warm absorber phenomenon, the primary goal 
of the work reported here is the outflow dynamics . 



2. Basic equations 

The equations to be solved include: 



|^ + V-(pv) = 0, (1) 

P + (V • V)V^ = -Vp - pV$ + p grad, (2) 

|! + V(v(e + p + p$)) = if, (3) 

these are the conservation equations for: mass, momentum and energy. Heating and cooling 
processes are described by the function if(erg cm"'^ s^^); e - is the sum of the kinetic and 
internal energy densities: e = pv'^/2 + e. Three components of the velocity are taken 
into account. However, the adopted azimuthal symmetry implies that all d/d(j) = 0. The 
gas is assumed to have a constant adiabatic index, so that e = P/(7 — 1) with 7 = 5/3. 
Equations iQ-^ are cast in a non - dimensional form with the characteristic scales set by 
the properties of the plasma orbiting at a distance of i?o = 1 PC from a black hole with mass 
Mg in units of IO^Mq: time in measured units of to(s) = 4.7 ■ 10^^ Vpi^ M^^^'^ , velocity in 
units of K)(cms-^) = 6.6 ■ 10^Mq^\~c^^. 



We anticipate that the conditions at the interface of a torus and an outflow driven by 
X-ray heating will resemble the X-ray self-excited wind in X-ray binaries. That is, an X-ray 
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heated skin is created at the outer layers of the companion star atmosphere as a result of 
X-ray heating by the compact object. The temperature there rises almost discontinuously 
to some value that allows the gas to escape the gravity of the star. The flow below the jump 
is subsonic and is determined by the details of the heating, cooling, and transfer processes in 
this region (IBasko et al. 19771) . In such a situation, the distribution of temperature depends 
on the heating and cooling, including both adiabatic effects and radiation, in the equations of 
hydrodynamics. This can only be calculated accurately if the interior of the star is included 
as part of the computational domain, and the same applies to the problem of interest to us: 
the torus X-ray excited wind. We, therefore include in our computational domain the entire 
volume containing the torus and we start our simulations from a rotating torus configuration 
in equilibrium in the external gravitational field of the black hole. 

A stationary, polytropic, toroidal structure with constant angular momentum was sug- 
gested by dPapaloizou fc Pringle 19841 ) (PP-torus), and was shown to be stable in 2.5 di- 



mensions against radial or axial perturbations, an found to be unstable to non-axisymmetric 
perturbations (this effect cannot be numerically investigated in axisymmetry) . In our dimen- 
sionless units the locus of the maximum torus density pmax occurs at Rq, so that pmax = 1- 
The distortion of the torus is controlled by the parameter d = {zu~^ + w~) / {2Rq), where 
refers to the inner or outer edge of the torus. No replenishing of the torus gas is assumed, 
i.e. gas can only leave computational domain. The drawback of this approach is that it is 
not possible to obtain a true stationary solution of the problem. Nonetheless, as we will 
show, the torus flow can be regarded as a quasi-equlibrium on timescales shorter than the 
depletion timescale, and can thus provide representative snapshots of the character of the 
flow. 



2.1. X-ray exited wind: evolution of the wind and the torus 

The thermodynamic properties of X-ray heated gas depend on the spectrum of the 
incident radiation as well as on the local atomic physics. Under the assumption of photo- 
ionization equilibrium the thermodynamic state of the gas can be parameterized in terms 
of ratio of radiation energy density to baryon density ([Tarter, Tucker, Salpeter 1969D: 



^ = 4 vr -Fx/n, where Fx = L^e ^/ (47rr^) is the local X-ray flux, Lx is the X-ray luminosity of 

the nuclei, t = K^p dr -is the optical depth. We assume that the attenuation is dominated 
Jo 

by Thomson scattering k, = 0.2(1 -|- Xh) — 0.4 cm^ g~^, where is the mass fraction of hy- 
drogen, and the factor - describes approximately the attenuation of the radiation on the 
way from the source toward a fiducial point. Assuming that the there is a fraction rj^ of the 
total luminosity Lbh available in X-rays and that the disk radiates a fraction rjedd of its Ed- 
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dington luminosity L^dd = 1-25 ■ 10^'' Mq we estimate: ^ ~ 4 ■ 10^ ■ rj^ rj^dd ^e/ (-^23 '"pc), where 
-is the distance in parcecs and A^23 is the column density in 10^^ cm~^. If the wind is driven 
by thermal evaporative mechanism, then its pressure and temperature are determined by the 
condition of photo-ionization equilibrium, if the dynamical time within the flow is much larger 
than the characteristic time of the heating-cooling. In a stationary flow the rate of the Comp- 
ton and photo-ionization heating should be balanced by the Compton, radiative recombina- 
tion, bremsstrahlung and line cooling respectively. In the calculations we adopt approximate 
formulas of ( IBlondin 199411 for these processes. The radiative acceleration, grad, accounts for 
both the acceleration due to lines (/rad = (Fuv/€/c) fct-" ((1 + - 1) /ri;,", where 

Luv is the UV flux, t = aeVthp/\dv/dr\, Vth is the thermal velocity, cxe is the electron- 
scattering opacity, and a = 0.5 ( [Castor, Abbott, fc Klein 1975 IStevens fc Kallman 1990p . 



and for radiation pressure on electrons Qcont = Fuyk/c. In our approach, the radiation 
force which is due to lines and also k and Tmax depend on the ionization parameter ^, and 
we adopt the description made by (IStevens Kallman 1990p . The ionization parameter ^ 
is calculated self-consistently at each grid zone, at each time step taking into account the 
attenuation of the X-ray flux due to flnite radial optical depth. 

We adopt spherical-polar (r, 6, 0) coordinates for calculation of a wind heated by a 
central source of radiation. We extend the computational domain from rin = 0.01 pc to 
'"out = 50 pc and the polar domain 9i from to vr with no assumption about the equatorial 
symmetry. Models presented here have been calculated with a N^x Ng = 100 x 100 grid. We 
have also performed test calculations on the 300 x 300 grid flnding no noticeable difference 
with the result presented here. We postpone the discussion of this and related issues to 
a separate publication. Both angle and radial grids are non-uniformly spaced: the angle 
gradually increases towards 6 = n/i with fewer angular points towards the pole (6* = 0) 
and equator {9 = n/2), and analogously spaced in the southern hemisphere. The use of a 
logarithmically spaced radial grid allows for the "throat" of the torus to be resolved. We 
set outflowing boundary conditions at both ends of the radial grid and axially symmetric at 
9 = 0, n. This initial torus is imbedded in a low density gas with p = I0~'^pmax, where pmax 
is the maximum density inside the initial torus. To solve the hydrodynamical part of the 
system of equations ([l])-([3]) we use the ZEUS-2D code described by ( IStone Norman 1992p . 



As a test we have evolved a torus distribution of matter for t = 2 (two rotational periods in 
our units) and found the whole conflguration stable. We include the heating term into the 
hydrodynamical equations using a fully implicit approach because the characteristic heating 
time scale can be much shorter than the dynamical time. 

We assume Mq = 1, rjedd = 0.4 , r/x = 0.5 and also that 0.5 Lbh is radiated in UV and 
contributes to the flow energetics through the radiation force grad- The initial torus conflg- 
uration chosen here is only mildly distorted, the distortion parameter: d = 1.25. Following 
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dPapaloizou fc Pringle 1984]) , we assume that the distribution of specific angular momentum 



inside the initial torus, is constant. To ensure that our initial torus is in equilibrium, we 



follow ( Papaloizou fc Pringle 1984 ) with the exception that we reduce the effective gravity 
GM-Bn/r'^il — fjcdd) due to the effect of the radiation pressure. However the real continuum 
radiation pressure, which is adopted in the calculations, is a factor of e"'^'' less because of the 
optical depth effect. Our initial torus is Thomson thick Tr{zu^,9 = 7r/2) ~ 4.35. Although 
the torus itself is not the primary topic of this study, the X-ray illumination does affect the 
interior and average properties of the torus. During the evolution the locus of the maximum 
torus density (located initially at ca = 1) shifts towards larger zu. This happens primarily 
because of the over-pressured region created by the evaporating gas and also because of 
radiation pressure. This allows for the energy to be deposited directly to the torus interior. 
In our case, direct radiative heating produces gradients of gas pressure in the torus interior. 
Thus, by making the initial torus marginally optically thick, we crudely mimic the model in 
which the torus is sustained against vertical collapse by the gradients of infrared radiation 
energy density (IKrolik 20070 . In what follows, we describe the torus outflow behavior at 
three different times after the initial turn-on of X-rays. 

After one orbital period, an over-pressured region extends to r ~ 4pc throughout the 
area that is not shadowed by the high density torus. At this time the non-shadowed area 
lies at 6' < 56°. In this space, an axisymmetric region exists between 0.7pc < zo < 2.5pc 
and z < 7 where the temperature is everywhere, of the order of the local virial temperature 
(except in the torus interior) : T ^ 5T^ij.{r). Closer to the black hole, the heating is much 
stronger, leading to T ~ 10''K. High temperature, low density gas fills the space in the 
"throat" of the torus. The outer edge of the torus has extends to ~ 5.3 pc in temperature, 
and to ~ 4.5 pc in density contours. The ionization parameter in the same region is ^ ~ 600. 
A wind is well-developed: the velocity field in the vicinity of the high density torus, follows 
the equal pressure contours; the maximum radial velocity is observed close to the axes at 
Vmax{(^ = 3°) = 1180 km s~^, the second local maximum is at 6* = 12°: fmax = 850kms~^. 
The maximum velocity gradually decreases towards higher inclinations fmax(^ = 27°) = 
343kms~^; and at = 60° it is fmax = 110 km s^^. Closer to black hole, at about 0.2 pc, the 
vg component of the velocity is more pronounced, until heights z ~ 1.4pc are reached forming 
a base of the evaporative wind. After this time, the flow is approximately symmetrical in 
both hemispheres. Because we are solving equations of ideal hydrodynamics (with only a 
small numerical viscosity), accretion through inner boundary (at r=0) is almost negligible: 
Min{MQ yr"^) < 10"^ The mass-loss rate M^^^(M0 yr'^ sterrad"^) peaks at ^ ~ 10°, i.e. at 
much higher inclinations than Wmax- The total mass-loss rate at this stage is M{Mq yr~^) ~ 
3- 10"^ 



Figure [T] shows the flow streamlines and density contours at three orbital periods. At 



- 7- 



this time the high pressure region expands vertically, to 2; ^ 10 pc from the inner torus. The 
torus inner edge can be inferred from the temperature and density maps: ~ 1.3 pc. Inside 
a throat of this torus, T ~ lOTvir. The contours of T become more flattened: 2;max — 4pc, 
^ 8 pc. The conservation of angular momentum helps to keep the outer edge of the torus 
at ~ 5pc. However, there is an envelope around this dense torus: 2;max — 4pc, ru^ ~ 8.1 pc. 
The values of the minimum ionization parameter and the column density vary significantly 
with the inclination angle: ^min(^ = 8°) = 2000, N22 = 0.3; at ^ ~ 25°, the ionization 
parameter decreases significantly: .^min(^ = 26°) = 63, N22 = 2.3; and at higher inclination, 
the ionization parameter is quite low: .^min(^ = 53°) = 6.6, and the obscuration is quite 
high: N22 = 75. The mass loss rate changes slowly dM / dt{{MQ yr~^)/yr) ~ 1.53 ■ 10^^. The 
mass-loss rate peaks at M^'^^ = 0.05 at 6* ^ 20°, where v^^^{e = 20°) = 616 km s'^ the 
maximum value of fmax(^ = 5°) = 1000 km s~^; The total mass-loss rate is M(M0 yr~^) ~ 
0.08. It is this evolved model that gives mass-loss rate in accord with the estimates of 
qKrohk fc Begelman 198"6t. 



Figure [2] shows the flow at five orbital periods. At this time the density maximum is at 
ru = 3.6. Remarkably, the inner edge of the inner, dense torus is at the same distance from 
the black hole: zu~ ~ 1.3 pc, as it was at three orbital periods, while the outer edge is at 
CtJ"*" ~ 6.2 pc. The position of the outer edge of the torus, as infered from the temperature 
maps, extends to 11 pc. Outside this region the flow temperature is T ~ 2.5Tvii.. The 
significant drop of ionization parameter from ~ 10'^ to ~ 10, occurs again a.t 6 = 25°, where 
the column also rises to a ~ 10^^. Aspect ratio of the torus: Ar/r ~ 1 in accord with what is 
inferred from observations ( Krolik fc Begelman 1986 : IJaffe et al. 2004|l . Almost everywhere 



in the wind the poloidal component of the velocity is determmined by Vr- However, closer to 
the black hole, the vg component is important, vg ~ Vr- The mass-loss rate is approximately 
constant at 15° < ^ < 60°: Mj^'^^ = 0.07, and has a peak at ^ = 74°: M^^^ = 0.16. The total 
mass-loss rate is M(M0 yr""*^) ~ 0.25 and c/M/dt ( (M© yr^^)/yr) ~ 10~^ . The maximum 
value of fmax(^ = 3°) = 1032 km s~^; another local peak is at = 15°, where the maximum 
velocity is fmax = 606kms~\ and the minimum Vma^{0 = 90°) = 180kms~^. The torus is 
losing mass in all directions, although with very different speed at different inclinations. 

We have performed multidimensional simulations of a wind excited by X-ray heating 
of a rotationally supported torus. Our dynamical calculations predict the velocity field and 
the spatial distribution of opacity and emissivity of the gas in both continuum and lines as a 
function of position. ( Synthesizing the spectrum provides a key test for our model ). We have 
calculated synthetic spectra by integrating over the wind and using the formal solution of the 
transfer equation. This is done by passing the results of the hydrodynamical calculations 
into a code which calls the subroutines from the xstar code flKallman fc Bautista 200"T1) 



to calculate the transmitted fiux. In this paper we present only absorption spectra. This 
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approximately represents the case when warm absorber gas lies on the line of sight to the 
observer and we see only transmitted spectra. These are shown in Figure [3} for various 
angles 9 = 26°, 35°, 44° at three orbital periods. These show that our simulations provide 
gas with intermediate ionization and sufficient column to produce observable warm absorber 
features. This supports the idea that it is this wind that is responsible for the warm absorber 
gas observed in Seyfert galaxies. More detailed results, including dependence of these results 
on different parameters as well as limitations of the adopted model and comparisons with 
observations, will be discussed in a separate paper. 

This research was supported by an appointment to the NASA Postdoctoral Program at 
the NASA Goddard Space Flight Center, administered by Oak Ridge Associated Universities 
through a contract with NASA, and by grants from the NASA Astrophysics Theory Program 
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LOG DENSITY AND VELOCITY AT 1=3 
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Fig. 1. — Velocity streamlines superimposed on contours of log(n/no) after three orbital 
periods. Note, that streamlines carry information about the direction only i.e. the magnitude 
of the velocity at low inclination is much higher than near the equatorial plane. 




Fig. 2. — Color-intensity plots of log(n/no) (left panel), and log(P/Po) (right panel) at five 
orbital periods; at northern hemisphere superimposed with velocity vectors. 



- 12 - 



9=26 
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Fig. 3. — The transmitted spectra observed at three orbital periods from 2.5D model 
of the evaporative torus wind. Each panel corresponds to a different viewing angle, 
6 — 26°, 35°, 44°; x-axis: energy; y-axis: normalized, transmitted flux. The initial spec- 
trum is a power law with energy index a — 1. 



